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We study analytically the plasmonic modes in the graphene-coated dielectric nanowire, based 
on the explicit form of nonlinear surface conductivity of graphene. The propagation constants of 
different plasmonic modes can be tuned by the input power at the order of a few tenths of mW. 
The lower and upper mode bifurcation branches are connected at the limitation value of the input 
power. Moreover, due to the nonlinearity of graphene, the dispersion curves of plasmonic modes at 
different input powers form an energy band, which is in sharp contrast with the single dispersion 
curve in the limit of zero input power. 


Nonlinear plasmonics is a newly developed but explosive growing field which not only offers extreme light ma¬ 
nipulation at the subwavelength scale [l|, but also provides an universal method to scale down the conventional 
nonlinear optical devices to the chip scale Q. Due to the strong local electromagnetic fields, nonlinear plasmonic 
effects can originate from the adjacent nonlinear dielectric media. In the past years, nonlinear plasmonic modes in 
metal-dielectric[^-@, dielectric-metal-dielectric [3-[IIl, and metal-dielectric-metal [H, [l^l planar structures have been 
studied extensively. Meanwhile, the existence of discrete solitons in nonlinear dielectric media embedded with periodic 
metallic films 0, nanowires [l^ . and nanorings [l^ have been proposed recently. 

As a newly discovered two dimensional electromagnetic material, graphene has received extensive attention in nonlin¬ 
ear plasmonics. Compared with ordinary dielectric media, graphene has a high nonlinear susceptibility [Tol - I^ . which 
is promising to release the demand of high input power in current nonlinear plasmonics. Until now, some basic phe¬ 
nomena based on graphene induced nonlinearities have been studied, e.g., solitons supported by monolayer graphene 
[ 2 ^ and multilayer graphene [25l - [?7j . However, little attention has been paid to the two dimensional graphene-based 
nonlinear plasmonic waveguides, although two dimensional structures are more favourable as fundamental building 
blocks of plasmonic waveguide arrays and plasmonic lattices [l^ . . 

In this Letter, we give explicitly the tensor form of nonlinear surface conductivity of graphene in the classical 
frequency range. As a simplest structure of the two dimensional waveguides, the plasmonic modes with different 
orders are presented analytically in the graphene-coated dielectric nanowire, where the nonlinearity of graphene is 
considered. Meanwhile, the dependencies between the propagation constants and the input power are discussed. 
Moreover, the energy band which is formed by dispersion cures at different input powers, is studied by taking the 
fundamental mode as an example. 

Due to the two dimensional nature of graphene, its third order surface conductivity is a fourth-order tensor with 
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16 elements. Previous studies mainly focus on the element ai^xx [H, HO, [13, whereas other elements are needed 
when the incident electric field has two components along the graphene surface. For this reason, we assume that 
the doped graphene monolayer is placed on the xy plane and a time-dependent electric field of the form E (t) = 
[Ex exp (—iujt) + C.C.] X -I- [Ey exp {—iojt) c.c.] ij is applied, where ui is the angular frequency of the electromagnetic 
field. Basicall y th e optical response of graphene is contributed by both the intraband and interband electronic 
transitions [221123. However, in the classical frequency range hut < /ic, namely the photon energy is smaller than the 
chemical potential, the optical response of graphene is dominated by the intraband transitions. Under the relaxation 
time approximation and neglecting the interband processes, the transport properties of electrons in graphene are 
governed by the following Boltzmann equation 




/(k,<)-/o(k) 


( 1 ) 


where k = {kx,ky) is the wave vector, /(k, t) is the nonequilibrium distribution function, /o (k) = 

1/ [l -I- is the equilibrium Fermi-Dirac distribution function, e (k) = vpJi^Jk'^ + ky is the Dirac cone 

spectrum of charge carriers in graphene, vp = c/300 is the Fermi velocity, c is the velocity of light in free space, 
T is the relaxation time, —e is the charge of an electron, h is the reduced Plank’s constant, kp is the Boltzmann’s 
constant, and T is the temperature. 

The exact solution of Eq. o at CUT ^ 1 is 


/(k,i) = 


-t/T 


dt'e* [k -f K (t, t ')], 


( 2 ) 


where k (t, t') = | E (t") dt”. The surface current along the graphene surface can be expressed as 

de (k) 


jW = -4 


(27r)^ h . 


dkf (k, t) 


dk 


( 3 ) 


where the factor 4 is due to spin degeneracy and valley degeneracy [27|| . In the low temperature limit T —>■ 0, the 
nonequilibrium distribution function can be replaced by the Heaviside step function. Thus the surface current reduces 
to 




ji (t) = - 


J —OO 


where 


h {t, t') = [ dk =H [yc- e{k + K {t, t '))], 

J ^Jkl + kl 


( 4 ) 


( 5 ) 


and i = x,y. Eq. ([S]) can be calculated by expanding the integrand with respect to Kx and Ky (up to the 
third order) and integrating over the Fermi surface. Calculation shows that 


li (t, t') = -kpTTKi {t, t') + -^Ki {t, t') K^ (t, t '), 

skp 


( 6 ) 


where kp = yc/vpTi is the Fermi wave vector. Without loss of generality, we consider the surface current in x direction. 
According to Eqs. (|H) and ([5]), and considering the equivalence between x and y coordinates, we obtain 


jx = jx (w) exp {-iujt) + jx (Sw) exp {-iSujt) + c.c. 


( 7 ) 


where the hrst term corresponds to the surface current with time dependence exp(—iwt) and the second term corre¬ 
sponds to the third harmonics. In what follows, we neglect the terms of third harmonics since the phase matching 
condition is required [3ll|. and obtain 


jxiuj) = 

+a^’::hEx\^ Ex. 


( 8 ) 

Thus crii = ay^J = 4(To^c/7’‘?i(l/T — iuj) is the linear surface conductivity, a^xyy = cr^yxy = cr^yyl = dy^yxx = cr^yxyl = 
Cyx’xy = (J^xxxl'i = cTy^’yy /3 = —^(Toe^v'jp/TTy,c^{l/T'^ + u}'^){1/t — Vluj) are the nonzero elements of nonlinear surface 
conductivity, and (Tq = j^Ti is the conductivity quantum. 
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FIG. 1: Structure of the graphene-coated dielectric nanowire. 


Under the limit of lot ^ 1, namely the relaxation time is large compared with the oscillation period of the incident 
electromagnetic field, the surface conductivity of graphene is 


where 


is the linear part of surface conductivity, 



(9) 

= 1 - 

Tlh, LO 

(10) 

. 

-1 - .2 ct 

(11) 

O'KflcTl 


is the nonlinear part of surface conductivity, and E|| is the electric field that is parallel to the graphene surface. Note 
in deriving Eqs. (|^- (fTT|l . we have used the condition E*Ex = E*Ey, which requires that the electric field E|| is 
linearly polarized. If the incident wave is circularly polarized or elliptically polarized, the tensor form of nonlinear 
surface conductivity should be used. 

Besides, Eqs. are valid under the approximation conditions of Tilo < fic and lot ^ 1. The carrier relaxation 

time T is determined by the carrier mobility /i as r = where the carrier mobility ^ of graphene film ranges 

from 1000 cm^/(V- s) in chemical vapor deposition (CVD)-grown graphene [2^ to 230 000 cm^/(V- s) in suspended 
exfoliated graphene [^. When using a moderate mobility of 10 000 cm^/(V- s) and with = 0.3 eV, the frequency 
of the electric field satisfies 0.5 THz <C lo/2'k < 45.3 THz. 

Considering the practical applications where two dimensional structures are more favourable, in the following we 
calculate the plasmonic modes of graphene-coated dielectric nanowire, where the structure is shown in Fig. [TJ For 
the modes propagating in the axial direction of the waveguide, the field can be expressed as 


A (r, 6, z,t) = A (r, 9) exp [i {Pz — cot)] + c.c.. 


( 12 ) 


where A denotes the electric field E or the magnetic field H, P is the propagation constant, and z is the propagation 
direction. From Eq. m, the z component of the field satisfies the following equation 


I d A, 1 d^Ax 

dj -2 r dr dO"^ 


(/ 3 ^ - k^e) A^ = 0 , 


(13) 


where ko = Lo^eofJ-o, and e is the relative permittivity of the material (e = ei for the inside dielectric nanowire and 
e = £2 = 1 for the outside air). Using the method of separating variables, the solutions for the m-th order plasmonic 
mode are 


E.{r,e)-- 

~ ^ -^m^m \ 

a ' 

^ a/ 

imO 

HAr,e)-- 

_ VPo/vo „ ^ 


) 

a 

V a> 


EAr,0) = 

j C K ( 

' r\ 

w — 

imO 

a ^ 

^ aJ 


Hx (r, 0) = 

. V^o/vo 

a 

1 w— 

V a 



(14) 

(15) 


(16) 


( 17 ) 
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for r > a, where a is the radius of the dielectric nanowire, u = a^/— kgSi, w = /3^ — kQS 2 , Am, Bm, Cm and Dm 

are the dimensionless undetermined constants, Pq = j ffg (E x H* + E* x H) • zdS is the input power, 770 = \/no/eo 
is the impendence of free space, and Im and Km are the m-th order modified Bessel function of the first kind and 
the second kind, respectively. Utilizing the above results, the other components of the electric field and magnetic 
field can be obtained from Maxwell equations [s^. Since the electric field that is parallel to the graphene surface is 
linearly polarized which satisfies E*Ex = E*Ey, we can use Eqs. ©-(HI]) to characterize the surface conductivity of 
graphene directly. Besides, since there are two components along the graphene surface, the contributions from both 
two components should been considered, which is different to the case of planar structures [IS, [IE 113 • According to 
the continuity conditions at r = a and at arbitrary 0, we can get the following equations 

Am^m (^) — CmKm ('^) , ( 1 ^) 

Bmljn (u) DjjiKfji {w') 

='i-^ [PmArnlm {u) + kouBmi'm (w)] i (19) 

— (w) -j- kQuBfjiI^ (^)] 

[l3mCmKm (w) + kowDmK'^ , ( 20 ) 

2 “1“ (^)] 

- \ [kQe2wCmK'^ {w) + PmDmKm (^^^)] 

= - t^Amlm (u) , (21) 


where 


cr = aA'i + |E|| {\Ez (a)l^ + \Eq (a)|^) , 

(o) = i Amlm (u) , 

a 

Eg (a) = i [PmAmlm (u) + kouBmi'm («)], 

and TO = 0,1, • • •. Meanwhile, since the input power is Pq; we obtain the normalization condition as follows 


( 22 ) 

(23) 

(24) 


l{e,Al + Bl)£ 

+ ‘^^AmBm {Slkl + [ Iml'mdr 

^ Jo 


rdr 


+ 


■Km. -^u Km.) rdr 


^ (aCj. + Dl) f 

Ja 

o poo 

■^—CmDm{e2kl+l3^) / KmK'^dr = l. 

^ J a 


(25) 


Thus, from the continuity conditions (fT51) - (IMl) and the normalization condition (1^51) . we can solve the five unknown 
parameters Am, Bm, Cm, Dm and /3 numerically. Specially, when the nonlinear surface conductivity of graphene is 
neglected, namely = 0, our nonlinear plasmonic waveguide reduces to the common linear plasmonic waveguide, 
and Eq s. (IT^ - (l25]) reduces to the dispersion relation of linear plasmonic modes, which have been discussed in Refs. 

iUllf. 

In what follows, we let / = a;/27r = 30 THz, /Xc = 0.3 eV, a = 100 nm, and ei = 3. Since the relaxation time of 
graphene ranges from 0.01 ps to 1 ps [IE, our parameters fulfill the approximation conditions of hcu < fie and wr ^ 1. 
As shown in Fig. [51 graphene-coated dielectric nanowire supports four orders of plasmonic modes, which correspond 
to TO = 0, 1, 2, and 3, respectively. Note when the input power tends to zero, the propagation constants reduce to 
that of the linear plasmonic modes [IE, where the nonlinear surface conductivity of graphene is neglected. In other 
words, the linear plasmonic modes in Refs. [sE are approximations of our nonlinear plasmonic modes in the limit of 
zero input power. 














5 



/>„(mW) 

FIG. 2: Mode bifurcation curves with different orders. The inset is the enlarged figure. The parameters are / = 30 THz, 
He = 0.3 eV, a — 100 nm, and ei — 3. 




FIG. 3: The distribution of plasmonic modes in the radial direction a.t 9 — 0 for (a) m = 0, (b) m = 1, (c) m = 2, and (d) 
m = 3, respectively. For comparison, the linear plasmonic modes are also plotted (black dashed curves). The parameters are 
/ = 30 THz, He = 0.3 eV, a — 100 nm, £i = 3, and Po — 0.2 mW. 


Starting from the zero power point, as the input power increases along the lower bifurcation branch, the field 
intensity at the graphene surface also increases. According to Eqs. (El-dlll), the surface conductivity of graphene 
decreases accordingly. Thus, the propagation constant of plasmonic mode increases with the enhancement of the field 
confinement, as shown in Fig. [3] For comparison, we also plot the corresponding linear plasmonic modes. Clearly, 
the propagation constants of these linear modes are equal to that of the nonlinear plasmonic modes in the limit of 
zero input power. 

In the lower bifurcation branch, the propagation constant increases monotonically with the input power, when the 
input power is below a limitation value. This limitation value is the maximum value of the allowed input power, 
where the mth mode only exists when the input power is below its corresponding limitation value. However, the 
field confinement can be enhanced further if the input power decreases from the limitation value along the upper 
bifurcation branch. Although the input power is decreased, the field intensity at the graphene surface is increased 
which insures the continuous growing of the propagation constant, as shown in Figs. [3]and[31 Note although there are 
intersection points when Pq = 0.339 mW for the curves with m = 0 and m = 1, and Pq = 0.213 mW for the curves 
with m = 0 and m = 2, degenerate states do not exist since these modes belong to different orders. Moreover, due to 
the high nonlinear surface conductivity of graphene, the propagation constants for different orders can be tuned by 
the input power at the order of a few tenths of mW, which cannot been realized by conventional nonlinear dielectric 
media [2^. Here we only show the four plasmonic modes which exhibit mode bifurcation at low input power, although 
other modes may also be supported by the graphene coated nanowire. 





























6 



FIG. 4: Dispersion relation of the fundamental plasmonic mode with the dependence of the input power. For comparison, the 
red dashed curve shows the dispersion relation of the corresponding linear plasmonic mode. The parameters are = 0.3 eV, 
a = 100 nm, ei = 3, and m = 0. 


Comparing the mode bifurcation curves with different orders, the limitation values for different plasmonic modes 
are not equal. At certain input powers, only the higher order modes exist, and the lower order modes vanish. This 
intriguing phenomenon is due to the self-action effect of graphene. Since the surface conductivity of graphene is 
dependent on the electric field that is parallel on the graphene surface, different field intensities are required to 
support the plasmonic modes with different orders, even at the same input power. Thus the nonlinear plasmonic 
waveguide may only support the higher order modes at certain input powers. 

For the linear plasmonic modes in Refs, [s^, the dispersion relation between the propagation constant and 
frequency / is a single curve, which is independent of the input power Pq. However, for the nonlinear plasmonic 
modes, the dispersion curves at different input powers form an energy band, as shown in Fig. |4l For simplicity, we 
only consider the fundamental mode with to = 0, where the incident frequency is tuned between / = 20 THz and 
/ = 45 THz to ensure the validity of the approximation conditions of hui < and wt ^ 1. 

Due to the nonlinearity of graphene, the plasmonic modes can only exist within a certain power range. For lower 
frequencies, the propagation constant can be tuned effectively by the input power. Whereas, the allowed band becomes 
narrow at high frequencies with the decrease of the limitation value of the input power. Actually, as the input power 
tends to zero, the dispersion relation reduces to that of the linear plasmonic mode, as shown by the red dashed curve 
in Fig. m Note the dispersion relation of the plasmonic mode would be different, if the realistic loss of graphene is 
considered. The further research is beyond the scope of this paper and we will show the relevant results else where. 

In conclusion, considering the vector nature of plasmonic modes in two dimensional waveguides, we derive the tensor 
form of nonlinear surface conductivity of graphene. The plasmonic modes with different orders are solved analytically 
in graphene-coated dielectric nanowire, where the propagation constant of each mode can be tuned effectively by the 
input power at the order of a few tenths of mW. The lower and upper mode bifurcation branches are connected at the 
limitation value of the input power. Moreover, due to the nonlinearity of graphene, the dispersion curves of plasmonic 
modes at different input powers form an energy band. Our work will provide important help to the research of other 
graphene-based nonlinear waveguides, especially plasmonic waveguide arrays and plasmonic lattices. 
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